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We show that when an isolated doublet of unbound states of a physical system becomes degenerate 
for some values of the control parameters of the system, the energy hypersurfaces representing the 
complex resonance energy eigenvalues as functions of the control parameters have an algebraic 
branch point of rank one in parameter space. Associated with this singularity in parameter space, 
the scattering matrix, St(E), and the Green's function, Gf + \k; r, r'), have one double pole in 
the unphysical sheet of the complex energy plane. We characterize the universal unfolding or 
deformation of a typical degeneracy point of two unbound states in parameter space by means of a 
universal 2-parameter family of functions which is contact equivalent to the pole position function 
of the isolated doublet of resonances at the exceptional point and includes all small perturbations 
of the degeneracy condition up to contact equivalence. 



PACS numbers: 03.65.Nk; 33.40.+I; 03.65.Ca; 03.65.Bz 



I. INTRODUCTION 

Recently, a great deal of attention has been given to 
the characterization of the singularities of the surfaces 
representing the complex resonance energy eigenvalues 
at a degeneracy of unbound states. This problem arises 
naturally in connection with the topological phase of un- 
bound states which was predicted by Hernandez, Mon- 
dragon and Jaureguip], and later and independently 
by W.D. Heiss@ and which was recently verified in a se- 
ries of beautiful experiments by P. von BrentanolT 
and the Darmstadt group |9|,ll0|, see also[n|. 



II. DEGENERACY OF RESONANCE ENERGY 
EIGENVALUES AS BRANCH POINTS IN 
PARAMETER SPACE 

In this short communication, we will consider the res- 
onance energy eigenvalues of a radial Schrodinger Hamil- 
tonian, Hr , with a potential V(r; x\ , x 2 ) which is a short 
ranged function of the radial distance, r, and depends on 
at least two external control parameters (x\,x 2 ). When 
the potential V(r;xi,x 2 ) has two regions of trapping, 
the physical system may have isolated doublets of res- 
onances which may become degenerate for some spe- 
cial values of the control parameters. For example, a 
double square barrier potential has isolated doublets of 
resonances which may become degenerate for some spe- 
cial values of the heights and widths of the barriers 

In the case under consideration, the regular and phys- 
ical solutions of the Hamiltonian arc functions of the ra- 
dial distance, r, the wave number, k, and the control 
parameters (xi,:^)- When necessary, we will stress this 
last functional dependence by adding the control param- 



eters (x\,x 2 ) to the other arguments after a semicolon. 
The energy eigenvalues £ n — (h 2 / '2m) fc^ of the Hamil- 



tonian H, 



are obtained from the zeroes of the Jost func- 
where k„ is such that 



tion, f(-k;x!,x 2 ) [!_, 

f(-k n ;xi,x 2 ) = 



(1) 



When k n lies in the fourth quadrant of the complex 
k— plane, Rek n > and Imk n < 0, the correspond- 
ing energy eigenvalue, £ n , is a complex resonance energy 
eigenvalue. 

The condition £[J defines, implicitly, the functions 
k n {xi,X2) as branches of a multivalued function 
which will be called the wave-number pole position func- 
tion. Each branch k n (xi,x 2 ) of the pole position func- 
tion is a continuous, single- valued function of the control 
parameters. When the physical system has an isolated 
doublet of resonances which become degenerate for some 
exceptional values of the external parameters, (x\,x 2 ), 
the corresponding two branches of the energy-pole posi- 
tion function, say £ n (xi,x 2 ) and £ n +i{xi, X2), are equal 
(cross or coincide) at that point. As will be shown below, 
at a degeneracy of resonances, the energy hypersurfaces 
representing the complex resonance energy eigenvalues as 
functions of the real control parameters have an algebraic 
branch point of square root type (rank one) in parameter 
space. 

Isolated doublet of resonances: Let us suppose 
that there is a finite bounded and connected region 
M. in parameter space and a finite domain T> in the 
fourth quadrant of the complex k— plane, such that, when 
(xi, X2) £ M., the Jost function has two and only two ze- 
roes, k n and k n +i, in the finite domain T> e C, all other 
zeroes of f(—k;Xi,x 2 ) lying outside V. Then, we say 
that the physical system has an isolated doublet of res- 
onances. To make this situation explicit, the two zeroes 
of /(— k; x±,x 2 ), corresponding to the isolated doublet of 
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resonances are explicitly factorized as 

f(-k;x x ,x 2 ) = (k - \(k n + k n+1 )) 2 - 
\{k n -k n+ i) 2 g n ,n+l(k,Xi,X 2 ). 



(2) 



When the physical system moves in parameter space 
from the ordinary point (xi,X2) to the exceptional 
point (x^yX^), the two simple zeroes, k n (xi,x 2 ) and 
k 7l +i(xi, x 2 ), coalesce into one double zero k d (x\, x%) in 
the fourth quadrant of the complex k— plane. 

If the external parameters take values in a neighbour- 
hood of the exceptional point {x\,x\) e M and k e T>, 
we may write 



9n,n+i(k; Xi,X 2 ) « g n ,n+i(kd, x*i,x* 2 ) ^ 0. 



(3) 



Then, 



1 1 1 

k - -^( k n(xi,x 2 ) + fc n+ i(a;i,X2)) 



xi,x 2 ) - k n+1 (xi,x 2 ) 



f(-k;xi,x 2 ) 

9n,n+l(kd] X* , X2) 



,(4) 



the coefficient [g 2)] multiplying 

f(—k;xi,x 2 ) may be understood as a finite, non- 
vanishing, constant scaling factor. 

The vanishing of the Jost function defines, implicitly, 
the pole position function fc„,n+i (xi, x 2 ) of the isolated 
doublet of resonances. Solving eq.10 for /c„ jn+ i, we get 

k n ,n+i{xi,x 2 ) = ^k n (xi,x 2 ) + k n+ i(x 1 ,x 2 ) S j + 

\J-^(kn{xi - x 2 ) - k n+ i{xi,x 2 )^ (5) 

with (x\,x 2 ) e M.- Since the argument of the square- 
root function is complex, it is necessary to specify the 
branch. Here and thereafter, the square root of any com- 
plex quantity F will be defined by 



F=\VF\exp(i~argF), 



< argF < 2ir (6) 



so that \VF\ = y/\F\ and the F— plane is cut along the 
real axis. 

Equation (JSJ relates the wave number-pole position 
function of the doublet of resonances to the wave number- 
pole position functions of the individual resonance states 
in the doublet. 

The analytical behaviour of the pole-position 
function at the exceptional point: 

The derivatives of the functions l/2(k n (xi, x 2 ) + 

k n +i{xi,x 2 )) and l/4(fc„(xi, x 2 ) - k n+ i(xi, x 2 )) 2 are fi- 
nite at the exceptional point. They may be computed 
from the Jost function with the help of the implicit func- 
tion theorem [l7| . 



d{k n {x- Li X%) — fc n +l(xi,X2)) 



dk 1 



k=k, 
df( — k;xi,x 2 ) 



(7) 



kd 



d [k n (zi,a 2 )+fcn+i (xi ,X2)J 
dxi 

/ X2 

f(~k;x 1 ,x 2 ) 
dxidk 



-1 



~ /(- ,& 2 ) 



{[( ! 



k=k d 



' f(-k;x 1 ,x 2 ) 
dfc 3 



k=k d 



k=k d 

7 d~f(-k;x 1 ,x 2 ) \ 1 1 
V dx ^ Jx 2 \k=k d i 



From these results, the first terms in a Taylor series ex- 
pansion of the functions l/2(fc n (xi, X2) + k n+ i(xi, x 2 )) 

and l/4(fc„(xi, x 2 ) — k n+ i{xi,x 2 )) 2 about the excep- 
tional point {x\,x 2 ), when substituted in eq.JSJ, give 

k n , n+1 (xi , x 2 ) = k d (xl,x 2 ) + Ak d (xi,x 2 )+ 

^\[c^\x 1 -xt) + ci 1) (x 2 -x* 2 )] (9) 

for (xi,x 2 ) in a neighbourhood of the exceptional point 
(x*, x 2 ). This result may readily be translated into a sim- 
ilar assertion for the resonance energy-pole position func- 
tion £ n , n +i{xi, x 2 ) and the energy eigenvalues, £ n {x\, x 2 ) 
and £ n+ i(xi, x 2 ), of the isolated doublet of resonances. 

Energy-pole position function: Let us take the 
square of both sides of eq.J^|, multiplying them by 
(h /2m) and recalling £ n = (h /2to)/c^, in the approxi- 
mation of ©, we get 



£ n ,n+i(xi,x 2 ) 



£ d {x* 1 ,x* 2 ) + A£ d (xi 1 x 2 ) 

(n,n+l{xi,X 2 ), (10) 



where 



e n , n+1 ( Xl ,x 2 ) = \/~[(R-£) + i(I-£)] 



(11) 



The components of the real fixed vectors R and / are 
the real and imaginary parts of the coefficients c| 
of (xi — x*) in the Taylor expansion of the function 

l/4(£ n (x 1 , x 2 ) — £ n+ i{x 1 , X2)) and the real vector £ is 
the position vector of the point {x\,x 2 ) relative to the 
exceptional point (x*,x 2 ) in parameter space. 
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= X1 
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-x* 2 ) 



R 



Re C\ 
Re C. 



(i) 



(i) 



Im C\ 
Ira C, 



(i) 



(i) 



(12) 
(13) 



The real and imaginary parts of the function 

(xi, x 2 ) are 



Ree 



n,n+l 



(xi,x 2 ) = ± 



1 



2^2 



(R-i) 2 +(f-i) 2 +R-i 

(14) 
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Ime nn+1 (x 1 ,x 2 ) = ± 



and 



2V2 



(15) 



1/2 



sign^Reen^+^jsign^Imen^+^j = signal '■ ttj (16) 

It follows from H14|) . that Ree n ^ n +i(xi, x 2 ) is a two 
branched function of (£1,^2) which may be repre- 
sented as a two-sheeted surface Sr, in a three di- 
mensional Euclidean space with cartesian coordinates 
(i2ee n) „ + i,£i,£ 2 )- The two branches of i?ee„ jn+ i(£i, £ 2 ) 
are represented by two sheets which are copies of the 
plane (^1,^2) cut along a line where the two branches of 
the function are joined smoothly. The cut is defined as 
the locus of the points where the argument of the square- 
root function in the right hand side of (|14l) vanishes. 

Therefore, the real part of the energy-pole position 
function, £ ntn+ i(x\, X2), as a function of the real parame- 
ters (xi,x 2 ), has an algebraic branchpoint of square root 
type ( rank one ) at the exceptional point with coordinates 
(x*,X2) in parameter space, and a branch cut along a 
line, Cji, that starts at the exceptional point and extends 
in the positive direction defined by the unit vector £ c sat- 
isfying. 



I-£ c = and R-Z C = -\R-£ C 



(17) 



A similar analysis shows that, the imaginary part of 
the energy-pole position function, Im £ n _ n+ i(xi,x 2 ), as 
a function of the real parameters (xi,x 2 ), also has an al- 
gebraic branchpoint of square root type (rank one) at the 
exceptional point with coordinates (x^a:^) in parameter 
space, and also has a branch cut along a line, Ci, that 
starts at the exceptional point and extends in the negative 

direction defined by the unit vector £ c satisfying eqs. \lT\) . 

The branch cut lines, Cr and Cj, are in orthogonal 
subspaces of a four dimensional Euclidean space with 
coordinates (i?ee„ jn+ i, Ime nn+ll £x, £2), but have one 
point in common, the exceptional point with coordinates 
(2;*, Xj). 

The individual resonance energy eigenvalues are con- 
ventionally asociated with the branches of the pole posi- 
tion function according to 



a 



= £((0,0) + A£„ in+1 (£i,&)- 

(w) 1 
R 2V2 



1(7 



(m) 1 



'I 2V2 
with m = n,n 

» 



\ .://-0 2 -(/-0- • i/<'-0 



1/2 

+ 
1/2 



(18) 



a 



1, and 



n+l 



» 



Re£ n — Re£ n+ i 
\Re£ n - Re£ n+ i\ 



Im£ n - Im£ n+1 
\Im£ n - Im£ n+ i\ 



(19) 



(20) 



Along the line £r, excluding the exceptional point 

(xl,x 2 ), 



but 



Re £ n (xi,x 2 ) = Re £ n+ i(xi,x 2 ) 



Im £ n (xi,x 2 ) ^ Im £ n+ i(xi,x 2 ). 



(21) 



(22) 



Similarly, along the line Ci, excluding the exceptional 
point, 



but 



Im £ n (xi,x 2 ) = Im £ n+ i(xi,x 2 ), 



Re £ n (xi,x 2 ) ^ Re £ n+ i(xi,x 2 ). 



(23) 



(24) 



Equality of the complex resonance energy eigenvalues 
(degeneracy of resonances) , £ n {x\,X2) = £ n +i(x*, x 2 ) = 
£ d (x\,x 2 ), occurs only at the exceptional point with co- 
ordinates [x*,x 2 ) m parameter space and only at that 
point. 

In consequence, in the complex energy plane, the cross- 
ing point of two simple resonance poles of the scattering 
matrix is an isolated point where the scattering matrix 
has one double resonance pole. 

Remark: In the general case, a variation of the vec- 
tor of parameters causes a perturbation of the energy 
eigenvalues. In the particular case of a double complex 
resonance energy eigenvalue £ d (x\, x 2 ), associated with a 
chain of length two of generalized Jordan-Gamow eigen- 
functions [l!|, we are considering here, the perturbation 
series expansion of the eigenvalues £ n ,£„+i about £ d in 
terms of the small parameter |£|, eas. <|18l20f) . takes the 
form of a Puiseux series 



£ n , n+l {x u x 2 ) = £ d (xt,x* 2 ) + \H\ 1/2 \/i[(R-i) + 
+A£ d (x 1 ,x 2 ) + o(\^ 2 ) 

(25) 

with fractional powers |£| J '/ 2 , j — 0, 1,2, ... of the small 
parameter |f| [Tllist. 



III. UNFOLDING OF THE DEGENERACY 
POINT 

Let us introduce a function fdoub(~ k; £1, £2) such that 
fdoub(-k;Zi,&) = \k- (fc d (0,0) + A«fc d (£i,6) 



and 



A^k d (x 1 ,x 2 ) = Y, 2 i= idPi 



(26) 



(27) 



Close to the exceptional point, the Jost function 
/(— k; £1, £2) and the family of functions fdoub(—k; £1, £2) 
are related by 

/(-*;;&, 6) « Ir-^fdoubi-k;^^) (28) 

9n,n+l\ K d, U, Uj 



4 



the term [g n ,n+i{kd, 0, 0)1 1 may be understood as 
non-vanishing scale factor. 

Hence, the two-parameters family of functions 



fdoub(-k;tu&) = [k- (k d + AMk d (£ 1) Z 2 ) 



(29) 

is contact equivalent to the Jost function /(— fc: £i , £2 ) at 
the exceptional point. It is also an unfolding 16, 20] of 
/( — fc;£i;£2) with the following features: 

1. It includes all possible small perturbations of the 
degeneracy conditions 



/(-fe;£i,£ 2 ) = 0, 



3/(-fc;£i,£ 2 ) 



dk 







/ g 2 /(-fc;£i,£ 2 ) 

V dk 2 
up to contact equivalence. 



^0 



(30) 



(31) 



2. It uses the minimum number of parameters, namely 
two, which is the codimension of the degeneracy Q- 
The parameters are (£i,£ 2 ). 

Therefore, fdoub{~ k; £1, £ 2 ) is a universal unfolding 
of the Jost function /(— fe;£i,£ 2 ) at the exceptional 
point where the degeneracy of unbound states occurs. 

The vanishing of fdoub{~ k; £1, £ 2 ) defines the approxi- 
mate wave number-pole position function 



fe„,n+i(£i,£ 2 ) = fc d +A^ ) „ +1 A: (i (£ 1 ,£ 2 )± 



r i(7?-£+,f-£) 



1/2 



(32) 

and the corresponding energy-pole position function 
£n,n+i(£i,£ 2 ) given in eq.JTUl)- 

Since the functions £ n (£i,£ 2 ) and £ n+ i(£i,£2) are ob- 
tained from the vanishing of the universal unfolding 
/dou&(-fc;£i,£ 2 ) °f the Jost function /(-fc;£i,£ 2 ) at the 
exceptional point, we are justified in saying that, the 
family of functions £ n (£i, £ 2 ) and f„+i(£i,£ 2 ), given in 
eas. flty) and I\iyi2u)) . is a universal unfolding or defor- 
mation of a generic degeneracy or crossing point of two 
unbound state energy eigenvalues, which is contact equiv- 
alent to the exact energy-pole position function of the iso- 
lated doublet of resonances at the exceptional point, and 
includes all small perturbations of the degeneracy condi- 
tions up to contact equivalence . 



IV. CROSSINGS AND ANTICROSSINGS OF 
RESONANCE ENERGIES AND WIDTHS 

Crossings or anticrossings of energies and widths are 
experimentally observed when the difference of complex 
energy eigenvalues £„(£i,£ 2 ) - £n+i(£i,£ 2 ) = AE - 
i(l/2)r is measured as function of one slowly varying 




0.08 



/n, s -0.04 



FIG. 1: The curves C n (iTi) and C n +i(jii) are the trajecto- 



ries traced by the points £n(£i, £2^) an d £n+i(£i, Q 1 ) 011 tnc 



t(ih 



hypersurface £n,n+i(£i, £2^) when the point (£i,£2 ±J ) moves 
along the straight line path m in parameter space. In the 
figure, the path tti runs parallel to the vertical axis and 
crosses the line Ci at a point (£i, c , £2 ) with £i, c < and 
^ < £2- The projections of Cni^i) and C n+ i(-Ki) on the 
plane (ImE,^) are sections of the surface Si; the projections 
of C n (ivi) and C n +i(iri) on the plane (Re£,£i) are sections 
of the surface Sr. The projections of C n (ji{) and C n +i(ji{) 
on the plane (Re£, Im£) are the trajectories of the S— matrix 
poles in the complex energy plane. In the figure, d — d* = £j 



parameter, £1, keeping the other constant, £2 = £2 . A 
crossing of energies occurs if the difference of real ener- 
gies vanishes, AE = 0, for some value £i jC of the varying 
parameter. An anticrossing of energies means that, for 
all values of the varying parameter, £1, the energies dif- 
fer, AE 7^ 0. Crossings and anticrossings of widths are 
similarly described. 

The experimentally determined dependence of the dif- 
ference of complex resonance energy eigenvalues on one 
control parameter, £1, while the other is kept constant, 



<£n(£l7£ 2 *' ) ) _ £n+l(£l>£ 2 ^) — £n,n+l(£l,£ 2 ^) 



(33) 



has a simple and straightforward geometrical inter- 
pretation, it is the intersection of the hypersurface 
£n,n+i(£i, £ 2 ) with the hyperplane defined by the con- 
dition 

To relate the geometrical properties of this intersection 
with the experimentally determined properties of cross- 
ings and anticrossings of energies and widths, let us con- 
sider a point (£1 , £2 ) in parameter space away from the 
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FIG. 2: The curves C n (n2) and C n +i{iT2) are the trajectories 
of the points £n(£,i,Q) and £n+i(£i, £f) on the hypersurface 
£n,n+i(£i> £2) when the point (^1,^2) nioves along a straight 
line path 7T2 that goes through the exceptional point (^1,^2) 
in parameter space. The projections of C n {~K2) and C n +i(n2) 
on the planes (Re£,(,i) and (Im£,( i i) are sections of the sur- 
faces Sr and 5/ respectively, and show a joint crossing of en- 
ergies and widths. The projections of C n {^2) and C n +i{^2) 
on the plane (Re£ , Im£) are two straight line trajectories of 
the 5— matrix poles crossing at 90° in the complex energy 
plane. At the crossing point, the two simple poles coalesce 
into one double pole of S(E). 



exceptional point. To this point corresponds the pair of 
non-degenerate resonance energy eigenvalues £ n (£1 , £2 ) 
and £„+i(£i, £2*''), represented by two points on the hy- 
persurface £„.„+! (£1, £2)- As the point (£i,£~2 ) moves on 
a straight line path 7Tj in parameter space, 



(34) 



the corresponding points, £n(£i,£ 2 ) ' ATL & £n+l(£i)£2 ) 
trace two curving trajectories, C n {iTi) and C'„+i(7ri) on 
the en,n+i(Ci; C2) hypersurface. Since £2 is kept con- 
stant at the fixed value £3 , the trajectories (sections) 
C n (iTi) and (7 n +i(7ri), may be represented as three- 
dimensional curves in a space £3 with cartesian coordi- 
nates (Ree, Ime, £1), see Figs. 1,2 and 3. The projec- 
tions of the curves C n (iTi) and C n +i(iTi) on the planes 
(i?ee,£i) and (7me,£i) are 



ite[6 TO (7ri)] = i?e£ m (£i,#) m = n,n + 1 (35) 



and 



Im[C m (wi)] = Im£m(€i, £2 ) m = n,n + 1 (36) 



FIG. 3: The curves C n (iT3) and CVi+i(7T3) are the trajecto- 
ries traced by the points £„ (£1 , £ 2 3 ' ) and Cj) 3 ') 011 the 
hypersurface £n.,n+i(£i, §2 3 ') when the point (£1,62 ) moves 
along a straight line path 7T3 going trough the point (£i, c , £ 2 3 ') 
with £i >c > £*. The path 713 crosses the line Cr. The pro- 
jections of C n {iT3) and C n +i{^3) on the plane (i?e£,£i) show 
a crossing, but the projections on the planes (Im£,£i) and 
(Re£ , Im£) do not cross. In the figure, £1 = d — d*. 



respectively. 

From eas. (|18l20[) . and keeping £ 2 = , we obtain 



A_B — E n — E n +i — (Re£ n — Re£ n +i 

<t<">V2 " 



( fl .0 2 + (/-0 2 + (^-0 



1/2-, 



and 



(37) 



Ar = ±(r„ - r„ + ij = Im[£ n +i) - (lm£ n J 



1/2 



^ (38) 

These expressions allow us to relate the terms (R ■ £) and 
(I ■ £) directly with observables of the isolated doublet of 
resonances. Taking the product of AEAT, and recalling 
ea. (fTf)|l . we get 



A£Ar = [I ■ £ 



(39) 



and taking the differences of the squares of the left hand 
sides of (l3*7Jl and we get 



S2— y 



(40) 
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At a crossing of energies AE vanishes, and at a cross- 
ing of widths Ar vanishes. Hence, the relation found in 
ea. lKfflJl means that a crossing of energies or widths can 
occur if and only if (I ■ 0c(«) vanishes 

For a vanishing (/ • £ c )jM = = AEAT, we find three 

cases, which are distinguished by the sign of (R ■ £ c )f<i> ■ 
From eqs. and 

1. (R ■ f c ) f(i) > implies AE ^ and AT = 0, i.e. 

^2 

energy anticrossing and width crossing. 

2. (R ■ |X (i) = implies AE = and Ar = 0, 

2 

that is, joint energy and width crossings, which is 
also degeneracy of the two complex resonance en- 
ergy eigenvalues. 

3. (R ■ |X W < implies AE = and AT ^ 0, i.e. 

^2 

energy crossing and width anticrossing. 

This rich physical scenario of crossings and anticrossings 
for the energies and widths of the complex resonance 
energy eigenvalues, extends a theorem of von Neumann 
and Wigner [2lJ for bound states to the case of unbound 
states. 

The general character of the crossing- anticrossing re- 
lations of the energies and widths of a mixing isolated 
doublet of resonances, discussed above, has been experi- 
mentally established by P. von Brentano and his collab- 
orators in a series of beautiful experiments Q, • 

A detailed account of these and other results will be 
published elsewhere [22II23]] 



V. SUMMARY AND CONCLUSIONS 



We developed the theory of the unfolding of the en- 
ergy eigenvalue surfaces close to a degeneracy point (ex- 
ceptional point) of two unbound states of a Hamiltonian 
depending on control parameters. From the knowledge of 
the Jost function, as function of the control parameters 
of the system, we derived a 2-parameter family of func- 
tions which is contact equivalent to the exact energy-pole 
position function at the exceptional point and includes 
all small perturbations of the degeneracy conditions. A 
simple and explicit, but very accurate, representation of 
the eigenenergy surfaces close to the exceptional point 
is obtained. In parameter space, the hypersurface rep- 
resenting the complex resonance energy eigenvalues has 
an algebraic branch point of rank one, and branch cuts 
in its real and imaginary parts extending in opposite di- 
rections in parameter space. The rich phenomenology of 
crossings and anticrossings of the energies and widths of 
the resonances of an isolated doublet of unbound states of 
a quantum system, observed when one control parameter 
is varied and the other is kept constant, is fully explained 
in terms of the local topology of the eigenenergy hyper- 
surface in the vecinity of the crossing point. 
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